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Abstract 

An entropic formulation of relativistic continuum mechanics is developed in the 
Landau-Lifshitz frame. We introduce two spatial scales, one being the small scale 
representing the linear size of each material particle and the other the large scale rep- 
resenting the linear size of a large system which consists of material particles and is 
to linearly regress to the equilibrium. We propose a local functional which is expected 
to represent the total entropy of the larger system and require the entropy functional 
to be maximized in the process of linear regression. We show that Onsager's original 
idea on linear regression can then be realized explicitly as current conservations with 
dissipative currents in the desired form. We demonstrate the effectiveness of this for- 
mulation by showing that one can treat a wide class of relativistic continuum materials, 
including standard relativistic viscous fluids and relativistic viscoelastic materials. 
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1. Introduction 

Thermodynamics is a macroscopic description of a physical system, taking the average of 
its microscopic degrees of freedom in both the spatial and the temporal directions. The 
spatial average is the usual coarse graining, while the temporal average is identified with 
a statistical (ensemble) average (see, e.g., [IJ). A global equilibrium of a given system is 
realized when the temporal average is taken sufficiently longer than any relaxation times of 
the system. This is characterized as the configuration which maximizes the total entropy. 
In usual observations, however, the time average is not taken for such a long time, so the 
materials we often encounter are not in global equilibrium and are regarded as being in the 
process of regression to equilibrium. 

We can say more about this regression to equilibrium if the microscopic interactions of 
a given system is sufficiently local in both the spatial and the temporal directions. Suppose 
that we divide the system into subsystems of small size (but still sufficiently large that 
thermodynamic descriptions work). Then, due to the locality, one may assume that each 
subsystem gets into equilibrium rapidly. This assumption is called the local equilibrium 
hypothesis. In the following, when writing spacetime coordinates x = (x^) we always assume 
that dynamical variables are already averaged around x in both the spatial and the temporal 
directions in such a way that this local equilibrium is realized well with required symmetries 
being respected. Each spatial unit over which the spatial average is taken is called a material 
particle. 

In the standard description of relativistic fluid mechanics, one adopts the formulation 
based on conserved currents [2JEI (for recent developments, see [HE] and references therein). 
For example, when describing a simple, relativistic fluid in a (D + l)-dimensional spacetime, 
one introduces as dynamical variables the energy- momentum tensor T^ u (x) , and the particle- 
number current (or the charge current) n^(x) ; in all, (D + 1)(D + 4)/2 degrees of freedom 
(= 14 when D — 3). They obey the following D + 2 conservation laws: 



Meanwhile, thermodynamic arguments based on the second law of thermodynamics (with 
a few natural assumptions) tell us that the (D + 1)(D + 4)/2 fundamental variables can be 
written only with D + 2 variables (and their spatial derivatives), and thus the conservations 
( II. ip are indeed enough for describing the dynamics of relativistic fluids. 

The method described in the previous paragraph is actually powerful and can be applied 
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to a wide class of continuum materials. However, it is rather difficult to see what actually 
happens entropically in the process of regression to the global minimum. In particular, the 
link between Onsager's idea of linear regression [61 [7] and the current conservations has not 
been given in a direct manner so far. (An attempt was first made by Casimir [8]. See also 

m 

In this paper, we propose a framework of linear nonequilibrium thermodynamics which 
directly realizes Onsager's idea of linear regression [6l [3, |8] by introducing the explicit form 
of the entropy functional which is local and to be maximized in the process of thermalization. 
We show that linear regression to the global equilibrium can be naturally represented in the 
form of current conservations with dissipative currents in the desired form. 

This paper is organized as follows. In Sec. 2, we give a relativistic theory of linear 
nonequilibrium thermodynamics, by closely following the idea of Onsager. We there intro- 
duce two spatial scales, e s and L s . The smaller scale e s represents the linear size of each 
material particle, while the larger scale L s represents that of a large system which consists 
of material particles and is to linearly regress to the equilibrium. We then propose the 
explicit form of the effective entropy functional and show that Onsager's original idea can 
be realized explicitly as current conservations with dissipative currents in the desired form. 
In Sec. 3, we demonstrate that usual relativistic fluid mechanics can be reproduced within 
our framework. In Sec. 4, we apply the formulation to a class of continuum materials [TO] . 
We propose a generally covariant generalization of the theory of viscoelasticity pUl HU [12] 
and show that the so-called rheology equations given in [TTJ [T2] can be naturally obtained 
in a generally covariant form. Section 5 is devoted to a conclusion and discussions. 

2. General theory 
2.1. Geometrical setup 

In order to define linear nonequilibrium thermodynamics in a generally covariant man- 
ner, we first make a few preparations. We consider continuum materials living in a (D + 
l)-dimensional Lorentzian manifold Ai . Its local coordinates are denoted by x M (// = 
0, 1, ■ ■ • , D) and the Lorentzian metric with signature (— , +, • • ■ , +) by g^ v {x) . We take the 
natural unit, h = c = k# = 1 , throughout this paper. In order to develop thermodynamics, 
we first need to specify a set of timeslices (or a foliation) on each of which we consider 
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spatial averages. In this article, we exclusively consider the foliation in which each timeslice 
is orthogonal to the energy flux [or the energy-momentum (D + l)-vector] p M , assuming 
that p M is hypersurface orthogonal. This choice of foliation is called the Landau-Lifshitz 
frame. We parametrize the hypersurfaces with a real parameter t as T, t [t G M.) and intro- 
duce a coordinate system x = = x) such that x° = t, x = (x l ) (i = 1 , • • • ,D). 
The unit normal to the hypersurfaces (called the velocity field) is given by 

u^x) = g^(x) p v (x)/e{x) = p^(x)/e(x) , (2.1) 

where e(x) is the density of the proper energy (rest mass plus internal energy), 



e(x) = yj-g^x) p„(x) p v (x) . (2.2) 

Here and hereafter, indices are lowered (or raised) always with g^ u (or with its inverse g^ u ). 
The induced metric on a D-dimensional hypersurface passing through x is expressed as 

h^x) = g pv [x) + u^x) u v (x) = g^x) + p ^ x ^ Pu ^ . (2.3) 

e [X) 

We define the extrinsic curvature K^ v of the hypersurface as half the Lie derivative of h pv 
with respect to the velocity field v? : 

K„ v = X - £ u h fiv = X - hfhj 7 (V p u a + V a u p ) , (2.4) 

which measures the rate of change in the induced metric h pu as material particles flow along 
u M . Note that this tensor is symmetric and orthogonal to u^, K^ v u v = . 



In the Arnowitt-Deser-Misner (ADM) parametrization [T3], the metric and t 
field are represented with the lapse N(x) and the shifts N l (x) (i = 1, ■ ■ • , D) as 



re velocity 







ds 2 = g^x) dx» dx v = -N 2 {x) dt 2 + h^x) (dx i - N^x) dt) (dx j - N j (x) dt) , (2.5) 
u = u»(x) d, = -— dt + — ^ di u,{x) da? = -N(x) dt) . (2.6) 

With a given foliation, we still have the symmetry of foliation preserving diffeomorphisms 
that give rise to transformations only among the points on each timeslice. Using this residual 



throughout the present paper, we write a contravariant vector field v^(x) (/i = 0, 1, ■ • • , D) in a con- 
cise form as v = v fJ '(x) <9 M , regarding d p simply as a basis of a vector space (i.e., the tangent space in a 
mathematical terminology). For example, Eq. (|2.6j) stands for the two equations, u° = 1/N and u % = N l /N 
(i = 1, • ■ ■ , D). We denote by u^V M the covariant derivative along the vector field v ~ <9 M , which now 
acts on tensor fields as a derivative operator. 
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gauge symmetry we can impose the synchronized gauge, N l (x) = , so that the background 
metric and the velocity field become 



ds z = g^(x) dx" dx u = -N 2 (x) dt z + h {j (x) dx l dx 1 , (2.7) 

(v) = (o • (K) = (o 3 • <2 - 8) 



u = u^x)d^ = T77 -^— = —, (2.9) 



d_ d_ 

N(x) di~!h 

where r is the local proper time defined by dr = N dt . The volume element on the hy- 
persurface through spacetime point x is given by the D-iorm \/hd D x = W det (hij) d D x . 
The volume element of the total (D + l)-dimensional manifold is given by \f—g d D+l x = 
Ny/hdtd D x. Note that even after taking the synchronized gauge, there remains a residual 
gauge symmetry of reparametrizing t that labels the timeslices. 

For generic coordinates (i.e., not necessarily the synchronized coordinates), we introduce 
the time derivative D/Di as the operation that satisfies the following conditions: 

(i) t = 1 , (ii) — u 11 = , (iii) Leibniz rule. (2.10) 

The first condition ensures that At (D /Dt) certainly measures the difference between a 
quantity on X t+ At and that on E t . The second condition ensures that if a system of small 
spatial region can be regarded as being static at time t, so can be the system at time t + At. 
An obvious solution is given by the Lie derivative along the velocity u = u^d^ multiplied 
with the lapse iV : 

^- = N£ U . (2.11) 

This also satisfies the condition (iii). In the following, the time derivative D/Dt will be often 
abbreviated as the dot. 



2.2. Thermodynamic variables 

As was discussed in Introduction, around each spacetime point x = (x° = t, x) on timeslice 
H t , we make spatial and temporal (or ensemble) averages over (D + 1) -dimensional regions 
whose linear sizes we denote by e s and e t , respectively (see Fig. [TJ. We assume that e s and 
e t are both much smaller than the curvature radius of the background metric g^ v , so that 
we can take a local inertial frame at each spacetime point x which can be regarded as being 
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Figure 1: Local thermodynamic equilibrium. 

flat at least within the extension of (e s , e t ) . This implies that thermodynamic properties 
of each material particle do not depend on the curvature of the metric when discussing 
thermodynamics in each material particle. Since the affine connection is not covariant, our 
assumption means that the covariant local entropy at x depends only on the local value of 
the metric, g^ v (x) . 

We assume that the local thermodynamic properties of the material particle at x (already 
in its local equilibrium) are specified by the set of local quantities (b A (x), c T (x), d p (x)) . Here 
c 1 (x) denote the densities of the existing additive conserved quantities C 1 . b A (x) denote 
the "intrinsic" intensive variables possessed by each material particle (such as strains), and 
d p (x) denote the rest "external" intensive variables which further need to be introduced to 
characterize each subsystem thermodynamically (such as the background electromagnetic 
or gravitational fields). In this paper, we distinguish density quantities from other intensive 
quantities, and construct, by multiplying them with the spatial volume element yh, new 
quantities which are spatial densities on each timeslice. For example, the entropy density s 
and the densities c 1 of conserved charges are density quantities, and for them we construct 
the following spatial densities: 

s = Vhs, c^Vhc 1 . (2.12) 

The local equilibrium hypothesis implies that the local entropy s(x) is already maximized 
at each spacetime point x and is given as a function of the above local variables; s(x) = 
s(b A (x), c 7 (x), d p (x)) . This functional relation is sometimes called the fundamental relation 
in the entropy representation. 

In the synchronized gauge, due to the relation d/dt = N(x)d/dr, the proper energy 
density e(x) measured with the proper time r is related to the energy density e(x) measured 
with time t as 

e(x) = N(x)e(x) . (2.13) 
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Note that e(x) includes the gravitational potential through the factor N(x) . Accordingly, 
the local temperature T measured with r (T = (ds / <9e) _1 ) is related to the temperature T 
measured with t (T = (ds / de)^ 1 ) through the following Tolman law: 

T(x) = N(x)T(x) . (2.14) 

2.3. Entropy functional and the current conservations 

A local thermodynamic equilibrium is realized only in each material particle of spatial size 
e s averaged for a period of time, e t . If we observe a material at spacetime scales larger than 
(e s , e t ) , we need to take into account the effect that the material particles communicate 
with each other by exchanging conserved quantities (e.g., the energy- momentum and the 
particle number). The second law of thermodynamics tells us that, if boundary effects 
can be neglected, this should proceed such that the total entropy of the larger region gets 
increased. 

In order to describe such dynamics mathematically, we first introduce another spacetime 
scale (L s , L t ) which is much larger than (e s , e t ) and assign to each spacetime point x = 
(x° = t, x) on the timeslice H t a spatial region S X [L S ] of linear size L s (see Fig. [2]). We then 




Figure 2: Time evolution of material particles in the large region S X [L S ] . 

consider the total entropy of the region S X [L S ]: 

S(t; E X [L B }) = I d D ys(b A (t,y), ?(t,y), d p (t,y)) . (2.15) 

«/ [Lb] 
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The irreversible motion of intrinsic variables a r (x) = (b A (x), c 7 (x)) at x will proceed toward 
the equilibrium state of the region T, x [L s ] . Due to the condition L s 3> e s , we can assume 
that the influence from the surroundings of the region T, x [L s ] is not relevant to the dynamics 
of a r (x) because x is well inside the region. An equilibrium state of the region E X [L S ] will be 
realized when the observation is made for a long period of time, L t , and can be characterized 
by the condition 

^y L -]) =0 (216) 

Note that the functional derivative is taken only with respect to a spatial, .D- dimensional 
unit in the functional. We denote the values of a r (x) at the equilibrium by a r Q (x\ L s ) = 
(b A (x; L s ), Cq(x; L s )) , and will call the procedure to obtain Oq(x; L s ) from a r (x) the dynam- 
ical block-spin transformation. One here should note that, since c^(t,y) are the densities of 
conserved quantities, the variations (12.161) with respect to c'-type variables should be taken 
with total charges kept fixed at prescribed values: 

I d^^^EC^E.tLJ). (2.17) 

A simple analysis using the Lagrange multipliers shows that the condition of global equilib- 
rium is expressed locally as 

ds 



A (x) = and h;(x)VMx)=0, (2.18) 
where j3j is the thermodynamic variable conjugate to c 1 that is defined by 

Pi(x) = ^(x). (2.19) 

We here make a few comments. First, the procedure to define the equilibrium values 
cio(x; £s) at point x is carried out over its own region E X [L S ] ; for a point y different from x , 
a^y; L s ) should be obtained by solving the equation 5S(t; Sj / [L s ])/5a r (|/) = 0. However, 
when y is well inside the region E x [L s ] , in the approximation of linear regression the value 
of a r Q (y\ L s ) can be regarded as the same with the value a^y; T, X [L S ]) that is obtained by 
solving the equation 5S(t; Y> X [L S }) / 5a r (y) = 0. 

arger than 
as well as 

that the boundary effects can be safely neglected. At the same time, we also take the 
2 Their fluctuations can be roughly estimated to be Aa5/(ag) ~ (e & /L s ) D / 2 . 



The second comment is about the spatial scale L s . We take L s sufficiently 
e s such that a^y; L s ) themselves can be treated as thermodynamic variables^ 
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spacetime scale (L s , L t ) not too much larger than (e s , e t ) such that the values of a r (y) 
= (t, y) G £ X [L B ]) do not differ significantly from those of a r (x; L s ) . The latter condition 
ensures that the fluctuation of a r (y) around a r (y; L s ) can be well approximated by the 
Gaussian distribution^ So long as L s is chosen in this way, the time evolution of the local 
variables a r (x) at time scale e t can be analyzed elaborately by decomposing a r (x) into the 
large-scale variable a r Q (x) and the small-scale variable (a — ao) r (x) , 

a r {x) = a r (x) + (a-a ) r {x). (2.20) 

In fact, the variables a r (x) evolve with scale (L s , L t ) , and their variations at time scale e t 
can actually be regarded as being negligibly small (see Fig. [3]) . In contrast, the variables 
(a — ao) r (x) represent the Gaussian fluctuations in the region E^fLg] , and their evolutions 
can be analyzed with the dynamics of spacetime scale (e s , e t ) in the linear approximation. 




We now recall Onsager's discussion on the linear regression of nonequilibrium system 
[3 E El [T] . We first decompose the time derivative of a r (x) , a r (x) = Da r (x)/Dt, into the 
reversible and irreversible parts: 

a r (x) = [a r (x)} rcv + [a r (x)} irr . (2.21) 

With the decomposition of a r (x) itself [Eq. f)2.20p ]. the reversible part can be written as 

[a r (x)} vcv = a r (x) + [(d - d ) r (x)] rev , (2.22) 



3 We also require that (L s , L t ) be much smaller than the typical scales over which the dynamics may 
change substantially (such as the typical size of a continuum material or the typical scale of the change in 
background fields). This ensures that the dependence of Oq(.t; L s ) on the scale L s is negligibly mild. 



S 



where the first term describes the reversible, isentropic motion of the large-scale variable 
Oq(x) , and the second represents the part of (d — cto) r (x) that does not affect the increase 
of the total entropy S(t; S X [L S ]) . Since in the current discussions a r (x) deviate from the 
equilibrium values a^x) only slightly, we can assume that the irreversible motion of the 
small-scale variables obeys the linear regression law, 

[a r (x)]- m = [(a - d o r(*)] irr = -!\ (a - ao)'(x) . (2.23) 

Here the matrix 7 = (7^) has only positive eigenvalues, and we have used the relation 
Oq(x) = [dg (a;)] rev , that is, [do(x)]i rr = 0. We then expand the total entropy S(t; H X [L B ]) 
around a r (x) to second order and write the difference between the total entropy for the 
configuration a r (x) , S(t; E X [L B ]) = J 2 < L , d D y s(b A (t, y), c 1 (t, y), d p (t, y)) , and that at the 
equilibrium a r (x) , S {t; E X [L S ]) = f Sx[Ls] d D y s(b$(t,y), c J (t,y), d p (t,y)) , as 

AS(t; E x [L a ]) = S(t; V x [L a ]) - S (t; V x [L a ]) 



\ [ d D y ( a -a ) r (3 rs (a-a y. (2.24) 



2 

Here (3 rs is a symmetric, positive semidefinite operator (which may include spatial derivative 
operators). The thermodynamic forces are then found to be 

^ 5_ASp ; S [L B ]) = _ _ 

and thus, substituting them into Eq. (I2.23p . we obtain the equations of linear regression: 

[a r (x)U = V s f.(x) = V s SAS ^J La]) (L rs = (7/3- 1 )") • (2.26) 

The so-called phenomenological coefficients L rs can be shown to satisfy Onsager's recip- 
rocal relation [61 U\ E] 

L rs = (_i)KI+KI L sr , (2.27) 
where the index \a r \ expresses how the variables transform under time reversal, a r (x) — > 



-l)l ar la r l 



x 







The Curie principle says that L rs can be block diagonalized with respect 
to the transformation properties of the indices (r, s) under spatial rotations and the parity 
transformation [9], that is, under the subgroup 0(D) of the local Lorentz group 0(D, 1) 



4 When the system is subject to external fields H which change the sign under time reversal (like magnetic 
field), the reciprocal relation is expressed as L rs (H) = (— l)' a ' l + l aS l L sr (—H) . 
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in local inertial frames. For example, when a r constitute a contravariant vector, (a r ) = 
(a 11 ) , the equations of linear regression should be set for each of the normal and tangential 
components to the timeslice through x : 

[K X )\L\( X ) 



[IV 



5 AS 

Sa u (x) 
5 AS 
8a v (x) 



(2.28) 
(2.29) 



where for a contravariant vector we define v± = {—u^u v ) v u and v 1 ^ = h^ u v u (and similarly 
for covariant vectors). Covariance and positivity further imposes the condition that and 
L|" / should be expressed as = Lj_ u^u u (Lj_ > 0) and = L\\ h^ v (L\\ > 0) , respectively. 

In Eq. (12.241) . the total entropy S(t] £ X [L S ]) is expanded around the equilibrium states 
a o(y) (y £ S X [L S ]) as a quadratic form in (a — ao) r (y) . However, due to the equilibrium 
condition (12.181) . it should be more efficient to express it as a quadratic form in {b — bo) A {i)) 
and ' V ' v f5i{y) , because we then can ignore the subtleness existing in the constrained 
variations with respect to c 7 -type variables. This consideration leads us to propose that the 
total entropy can be effectively expressed by the following local functional (to be called the 
entropy functional): 



S(t; E X [L S }) 
= So(t; Tl x [L 8 \) 



A D yN^Ub-b li ) A V„ft 




it 



u, 







1 fil, vj 



U 



.) 



(2.30) 



Here the coefficient (^^fi j is a symmetric, positive semidefinite matrix. We have 
extracted the factor N from the coefficient matrices for later convenience (recall that N\fh = 
y/—g). The matrix elements in principle could be calculated from the fundamental relation, 
s = s(b A , c 1 , d p ) , once the foliation is fixed [see the Appendix |A] for the derivation of (12.301) 
for simple cases]. However, as we see below, their explicit forms need not be specified for 
the following discussions. 

With the entropy functional (I2.30p . the thermodynamic forces in the equations of linear 
regression Eq. ( 12.261) are calculated as 

5AS 



5b A 



~9 [£ab (b - b ) B + tJi V M ft - s M V M « (b - b ) B + l^ vJ V u f3j)} 
~g[^AB{b-b G ) B +^'V / A], (2.31) 
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fi = 6 -^ = V=g (&5/d&ds T ) v M [if (b - h) B + tr 3 "* v„p K ] 

~ Ns^j V„ [if (b - b ) B + i^ K V V /3 K ] , (2.32) 

where sai = d 2 s/db A dc I , s/j = \/hd 2 s/dc I dc J and s°j is the value of s/j at the equilibrium. 
In obtaining the last line of Eq. (I2.3ip we have neglected higher-order terms in the derivative 
expansion. Then the irreversible evolution of (b A , c 1 ) is given by 

[b\l = ^L AB f B + L AI (2.33) 
[5\l r = L IA f A + VhL IJ fj, (2.34) 

where we have multiplied (\fh) ±1 to make L rs tensors (not tensor densities). By further 
assuming that L IA = L AI = and by using Eqs. (12.311) and f 1 2 . 3 2 [) . the above equations can 
be rewritten as 

b A - [b] f ev = [6] i = -N L AB [i BC (b-b ) c + Fj ] , (2.35) 
c 7 - [£]L= ML = v^L IJ s° JK V»[if (b-b ) B + P K '» L V u p L ] . (2.36) 

The coefficient L IJ s° JK can be regarded as being constant in our derivative expansion, and, 
as we see in the next section, for the cases of interest the quantities c / = N£ u (\/hc I ) can 
be always written as y/—g V /i (c / w At ) . We thus can rewrite Eq. ( I2.36P in the form of current 
conservations: 

V M J /A1 = . (2.37) 

Here the current J Ifl consists of three parts: 

J*» = C V + jg + . (2.38) 

The first is the convective part, and the third is the dissipative part defined by 

= -L IJ s° JK [if (b - b ) B + t K ^ L V v p L ] • (2.39) 

The second comes from the equation [c]^ ev = —yf—g V M j!^ , which should be derived from 
another theory describing reversible, isentropic processes. For example, for simple viscous 
fluids, such theory should be that of ideal fluids. 

When there are only c J -type variables which are all scalar functions, the coefficients 
fyi, v j s hould take the form i^' uJ = -ih^ v {sq 1 ) ij (£ > 0) due to the positivity and the 
covariance, and the dissipative current is expressed as 

= £ L IJ W v V u /3j = i L IJ hT V„ (|^) ■ (2.40) 
11 



The second law of thermodynamics says that S(t; T, X [L S ]) is a monotonically increasing 
function of time t . Indeed, the entropy production rate P(t; H X [L S ]) in the region S X [L S ] is 
always positive for time scale e t due to the positivity of the phenomenological coefficients: 

P(t; Z X [L S }) = ^-AS(t; E X [L B }) = [ d D y 5 -^- if 

= A77 a W = / d D y L rs — — — — > . (2.41) 

Here we have neglected the contributions from d p because they would be of higher orders. We 
also have used the identity J E , L , d D y(SAS/Sa r ) [a r ] rcv = . One can further show that the 
entropy production is a monotonically decreasing function when [d r ] rev = 0; P(t; H X [L S ]) < 
0. 



3. Relativistic fluid mechanics 

3.1. Relativistic fluids in the Landau-Lifshitz frame 

We define relativistic fluids as continuum materials whose thermodynamic properties can be 
characterized only by the local energy-momentum p a = Vhp a (a — 0, 1, • • ■ , D) , the local 
particle number h = y/h n , and the background metric g^ u . 

To relativistic fluids, the formalism in the previous section can be applied with the 
following identifications: 



b A 


(ds/db A ) 


c 1 


Pi 


d p 


{d~s/dd p ) 


nothing 


nothing 


Pa 

h 


-u a /T 
-fx/T 







We here explain the entities in the list. We shall assume that s depends on p a only through 
the local energy e(p a ,g^ u ) = \f—g^ p^ p u so that the fundamental relation is expressed in 
the following form: 

s(p a A,9^u) = v(e(p a ,9r,),n,gp,) ■ (3.1) 

Then, by introducing the temperature T , the chemical potential /1 and the quasiconservative 
stress as 

dd{e,h,g,J) _ 1_ da(e,n,g^) _ // da(e,h,g^) _ Vh ^ 

de ~ T ' dn ~ T ' d 9flv ~ 2T ^ ' 1 } 
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and by using the identities 

dt 



P 



de 



9p a 



— u 



dg 



fJ,V 



~2<r 



eu^u 



(3.3) 



the variation of the local entropy s is expressed as 

5S = 



— S Pa - - 5n + — Tgj , 



where T/\ is the quasiconservative energy-momentum tensor: 



rpflU 

J (q) 



p^u u + r,^ 



(q) • 



(3.4) 



(3.5) 



In order to make this decomposition unique, we require that rf\ be orthogonal to : 



v - 

(q) 17 ~~ ' 



(3.6) 



Simple fluids (that have no specific spatial directions) are realized by taking r/^ = P h^ u 
with P the pressure. Then, by using the identity 5g iiv = h? v Sh^ = 25y/h/Vh, Eq. (EOI) 
becomes the standard expression for the local entropy s = a(e, h, g^ v ) of simple fluids, 5s = 
5a = (1/T) 5e — {n/T) 6n + P 8\fh] or equivalently, 5s = 5{s/ \fh) = (1/T) 5e — (p/T) 5n . 

As for the entropy functional, from the general expression (I2.30p without 6^-type vari- 
ables, we have 



AS 



d D yNVh 



ds 



ds \ 



+ m^dj^)d u 



dn) 



ds 
dn 



(3.7) 



We assume that the coefficient functions v g and m} iv are all orthogonal to the velocity 
field; i^a^p u v = , u 13 = , and m^ u u v = , and also that the second derivatives of the 
form d 2 s/dp a dn are small and can be neglected. The maximum of the entropy functional 
is then given by (p°,n°) that satisfy the following equations: 



ds \ 

.dp /3 

-MS) 



u 



~T 



n 
0. 



hT h a0 V v tf 







(3.8) 
(3.9) 



Note that these equations hold only within the spatial region S x [L s ] . 
The variation of the entropy functional (13. 7p is given by 



5AS 
5px 



- d 2 s 
9 dp x dp a A 



_ d 2 s 
9 dp x dp a A 



(3.10) 



13 



5AS 

5n 



d 2 s 



dn 



~2 v m 



d 2 S 



dn 



~2 v ^ 



m? v d v - 



77 J 



(3-11] 



By using the decomposition of the matrix d s/dp^dpy (negative-definite for each irreducible 
component) as 



d 2 s 



dp^dpv 



-c± v?u v - ci, hT 



(cx = -VSg(>0), q, = ^(>0)) 



(3.12) 



the vector 5AS/5p\ is decomposed as 

<5AS 



5AS 
5p\ 



+ 



5AS 

5pA 



with 



5AS' 
"5AS 



-u x u u ) ^ = iVc ± (-ttV) V 



8pp 



i' 



= h 



5 AS 

Spu 



New h Xu V 



rp C V, cr V P a 



(3.13) 

(3.14) 
(3.15) 



With the variations of the entropy functional Eqs. ( 13. lip . ( I3.14p . and A3. 15j) . we set the 
following equations of linear regression: 

~5AS~ 



[Pa ~ [Pa}v CV ] ± = [P. 
\_Pa [Pctjrev] m 



n — \n\ 



Pa\h 



n 



VhLj_ u a u\ 

\ThL\\ h aX 
5AS 



5px 
SAS 
Spx 



VhM 



5n 



(3.16) 
(3.17) 
(3.18) 



Note that the time derivatives of p a and n are given by 

^p a = N£ u p a = V^gV^(p a u^), (3.19) 

^-f l = N£ u n = ^\/, l (nu fl ). (3.20) 

Equations ( 13 . 1 9[) and ( I3.20p can be shown by using the identities £ u \/h = y/hW^u^ and 
P m Vq,m m = 0. The former identity can be proved as £ u \fh = (l/2)\/~hW£ u h fll/ = (1/2) 
Vhh^£ u9fMU = {lfflVhh^iVnUv + VvuJ = (l/2)Vhg^(V^u u + V v u M ) = ^V/. 
The latter follows from the facts that p [1 = eu^ and u^u^ = — 1 . One can easily show that 
the left-hand side of Eqs. (13. 16[) — ( T3 . 1 8[) can be rewritten as 
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(3.21) 



[Pa ~ [Pa]rcv]|| = h£ [y/^g 'V ^{Vv « M ) ~ \pu)rev] 



n — \n\ 



(3.22) 
(3.23) 



Furthermore, one expects that the reversible (or isentropic) evolutions of and n can be 
identified as the evolutions for ideal fluids, 



[Pc 



(3.24) 



The isentropic current Ta^ equals the quasiconservative stress tensor M , as is shown in 
the next section. Then Eqs. fl3.16p - fl3.18p become 



-u a u v ) V M T^ = - (-u a u v ) c ± L ± V, 



h:v fl T^ = -h:c ll L ll v, 

V M (nu M ) = Vh-^MV^ 



ft 1 p v u a 

rp c V, a v p a 



(3.25) 
(3.26) 
(3.27) 



We see that general covariance requires that c±L± = cy Lu (= Lp) . By assuming that Lp 
and L h = \fh (—d 2 s/dh 2 ) M can be regarded as being constant at this order of derivative 
expansion, the linear regression is expressed as current conservations: 



V M T/ = 0, V/ = 0, 



(3.28) 



where the energy-momentum tensor and the particle-number current are given by 

T/ = T^ )p + 4 d) ^ = e u v u p + + , n p = nu p + u p (3.29) 



with 



( d )M = L-F p V u u 

— rp C V, cr V P U 



L - m n»dj-E 



(3.30) 



For simple fluids that have no specific directions, the quantities given above can be 
parametrized as 



Tgj = P h pv (P: pressure) , 



- L p F u ' pa = C hT h pa + T) [h w h va + hT h up - {2 /D) hT h pa ] 



T 



L* m pu = aT 2 h pv 



(3.31) 
(3.32) 
(3.33) 



where ( , rj , and a (> 0) are, respectively, the bulk viscosity, the shear viscosity, and 
the diffusion constant. Note that the contributions from the rotation h p h"V[ p u a ^ have 
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been discarded since it vanishes in the assumption that the velocity field u p is hypersurface 
orthogonal. Then Eqs. fl3.28p - fl3.30p become the well-known conservation laws for viscous 
simple fluids in the Landau-Lifshitz frame: 

V M T^ = , VX = V M (n u» + v») = , (3.34) 
= e u»u v +{P- CV p u p ) W v - 2r]V {lM u u) , (3.35) 
/ = ^r^(^). (3.36) 

Here and hereafter, we define = (1/2) h p p h v a (A pa + A up - (2/D) h aP A aP h pa ) for a 

tensor A^ v . 



3.2. More on the entropy production 



The entropy production rate can also be calculated in the following way once we know that 
the current conservations f l3.34p - fl3.36p hold. First we consider the increase of the entropy 
functional during the time interval At (> e s ) : 

"t+At 



S(t + At; Z X [L S }) - S(t; E X [L S }) = J dtS(t; Z X [L S }) 
d D ys(p ll (t,y), h(t,y), g^{t,y)) 



t+At 

dt 

t J^ x [L Sj 
t+At p 

dt / d D y 

t+At 

d D+1 y 



ds i ds x. ds . 
k^Vv + -^r n + - — g^ 



(3.37) 



Since V v T pu = for T pu = p p u u + = p p u v + r/Jj + rgj , V M n^ = for n p = nu p + i/" 



and V u u^ = V v u^ 



Ufi ^ 7 ^ r (q) > the above equation can be rewritten as 



S(t + At; E X [L S ])-S(t; E X [L B }) 

"t+At 

d D+1 yv~9 



JT, X [L S ] 
t+At 



lD+1 



yv-g 



-u p 
. ~T 


MV- 


> 






T 



(3.38) 



The first equality means that t u u equals v in the absence of dissipation, and thus shows 



that the equality rp, 
further rewritten as 



holds for fluids. By integrating by parts, the equation can be 



5(t + At; 
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t+At 



i 



rp ' (d) 



T 



T 



(3.39) 



The last term on the right-hand side is the effect from the surroundings so that the entropy 
production in the bulk is given by 



S(t + At; E X [L B ))-S(t; £,[L S ]) 



l-D+l 



t+At 



yv-g 



rp T (d) ^ U U H 



£ 

T 



T 



(3.40) 



Since s can be expressed as / = gV /1 (su' 1 ) , the left-hand side of Eq. (I3.40p can be written 
as ft +At J E y L j d D+1 y y/— g V M [s u M — (fi/T) z/ m ] . Thus, by defining the entropy current as 
s M = — (/i/T) z/ M , Eq. (I3.40p can be expressed in a local form as 
1 



~T 



ji (d) "/w 1 " 

T 



£ 
T 



(3.41) 



Equation (13.41 \ might seem to be inconsistent with the entropy production rate given 
in Eq. ( I2.4ip . because Eq. (I2.4ip predicts that the entropy production rate is the spatial 
integral of a quadratic form in V^i^ and V M r^ . However, they are actually consistent. In 
fact, suppose that we expand u^/T and \xjT in Eq. (I3.38P around their equilibrium values 
(u%/T )(y; L s ) and (fi /T )(y; L s ) as 



£ = £) 
T T 



£ _ £) 

T T 



(3.42) 



Since the deviations /T — Uq/Tq and fi/T — hq/Tq are proportional to their thermodynamic 
forces, SAS/Spfj, and 8AS/8n, respectively, they can be expressed in the following form: 



T T ~ V " 



(3.43) 



T T 

where ^^ u ' pa and 7 A " / are some positive semidefinite matrices. Substituting them into Eq. 
( I3.38p . we have the following equation: 



Jt, x \l s ] w o Jo 

-V„[(r»>r°/T)K p(T ) V A rW A - V M [-r^(A*/r)] V, 



LZvs 



(3.44) 



Here (uq/Tq) t^ u can be neglected at this order of approximation because «{{r ( 



M ^(d) i/ 



U 1 * T, 



0, while Vjy[(/^o/T ) z/ 17 ] can be replaced by V u [(fi/T) u u ] . Thus, the second 
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term on the left-hand side of Eq. (I3.44p can be written as — J s , L , d D y yf—gV u [(fi/T) u u ~\ , 
so that the left-hand side of Eq. (13.441) agrees with that of Eq. (I3.40p . We thus find that the 
entropy production rate in the bulk can also be written as a quadratic form in V M ^ M and 



4. Relativistic viscoelastic fluid mechanics 



4.1. Intrinsic metric 

Once the motion of a material is specified by a velocity field u^{x) , the shape of the timeslice 
normal to is represented by the induced metric (12.31) . = + u^u u . If the relaxation 
time for the deformation of shape can be regarded as being very short, then, as discussed 
in Sec. [3j the fundamental relation, s(x) = a(^e(x),h(x), g fiu (x)) , does not depend on the 
hysteresis of the shape h^x) . Almost all materials (besides fluids), however, do not have 
such short relaxation times for the deformation of shape, and the time evolution of shape 
may not be described as a Markovian process; all the preceding history needs to be known 
in order to predict the future behavior of the system at a given initial timeslice. One then 
needs to extend the formalism such that the equations contain time derivatives of higher 
orders. This is out of the scope of the standard nonequilibrium thermodynamics. 

However, a class of materials still allow a standard thermodynamic description even when 
their relaxation time of shape is not short. This is performed by replacing the introduction 
of higher time derivatives with that of extra dynamical degrees of freedom. Viscoelastic 
materials considered below belong to this class of materials. 

According to the definition of Maxwell, viscoelastic materials behave as elastic solids at 
short time scales and as viscous fluids at long time scales (see, e.g., Sec. 36 in [H]). In order 
to understand how such materials evolve in time, we consider a material consisting of many 
molecules bonding each other and assume that the molecules first stay at their equilibrium 
positions in the absence of strains (as in the leftmost illustration of Fig. H]) [TTJ [12]. We 
now suppose that an external force is applied to deform the material. An internal strain 
is then produced in the body, and according to the definition, the accompanied internal 
stress can be treated as an elastic force, at least during short intervals of time. However, 
if we keep the deformation much longer than the relaxation time (characteristic to each 
material), then the bonding structure changes to maximize the entropy, and the internal 
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deformation 




Figure 4: Processes of deformation and strain relaxation. 

strain vanishes eventually as in the rightmost of Fig. |H The point is that two figures (the 
central and the rightmost) have the same shape (same induced metric) h^ v , but different 
bonding structures. 

In order to describe the internal bonding structure, we introduce at each moment another 
dynamical variable to be called the intrinsic metric^ which measures the shape that the 
material would take when all the internal strains are removed virtually [TQl [HI EE2]- To 
define such states, around each spacetime point x , we consider a small spatial subregion 
whose linear size A is much smaller than L s and enclose the subregion with an adiabatic 
and impermeable wall of vanishing tension (see Fig. [5]). We then cut it out of the bulk 




Figure 5: The procedure to obtain the intrinsic metric. 

S X [L S ] and leave it for a long time (comparable to L t ). The subregion will then undergo a 
relaxation to reach its equilibrium state. This causes a diffeomorphism in (D+l)-dimensional 
spacetime from the subregion to another subregion, mapping the original (or real) positions 
of material particles to their new positions in the virtual equilibrium, fx : x i— >■ f\{x) . Since 
the absolute value of the extrinsic curvature gets reduced in the course of relaxation, this 
diffeomorphism does not necessarily preserve the foliation, and the velocity field u^(x) will 
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change as u^(x) — > u'^x) .We then consider the pullbacks of the resulting fields in the virtual 
equilibrium: 

u^x; A) = (fxu')v(x) , g^(x; A) = (fx*g)^(x) , (4.1) 

and define the intrinsic velocity field u^x) and the intrinsic metric g^ v {x) as the values in 
the limit A — > : 

M x ) = J™ (/>%W , 9»»{x) = Jim (fxg)»»(x) ■ (4.2) 

Since g^ u u v is no longer orthogonal to the original hypersurface, the intrinsic metric 
g^vix) may have tilted components. We parametrize g M „ as 

g^v = - (1 + 20) u^Uy - e^u v - e u u^ + hfj, v 

(e/ = 0, h ia ,u v = 0, S F «" = 0). (4.3) 

The quantity vl + 20 represents the ratio of the temperature in the presence of strains 
to that in the absence of strains, because the Tolman law (T = NT = NT) tells that 
T/f = N/N = VI + 20, that is, = (T 2 — f 2 )/2T 2 ~ (T — f)/T . The quantity e„ 
represents the relative velocity of the real path of a material particle to that of its virtual 
path. The spatial metric h^ v is the same with the "strain metric" introduced by Eckart to 
embody "the principle of relaxability-in-the-small" in anelasticity pU]. (This was reinvented 
in [11] in the light of the covariant description of viscoelasticity under the foliation preserving 
diffeomorphisms.) We further introduce the strain tensor 

E^u(x) = - (g^ u {x) - g^ v {x)) 

= 9 u^u v + - (e M u„ + e v it J + , (4.4) 

where 

£nu{x) = g (V^) ~ h^ v (x)) (4.5) 

is the spatial strain tensor. The explicit forms of h^ u and h^ u under various deformations 
can be found in [12]. Note that if we define the extrinsic curvature associated with the 
spatial intrinsic metric h^ v as 

K,,,, = - A V = - (u x d x h^ u + h Xfl + d u u x h^x) , (4.6) 
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the following identity holds: 

£ e = K — K 

We denote the contraction of a spatial tensor A^ u with g^ u by tr A so that 
tr e = g^e^ = h^e^ , tr K = g^K^ = hTK^ . 



(4.7) 



(4.8) 



4.2. Linear regression and the rheology equations 

We now apply the framework of section [2] with the following identifications: 
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Pn 


-u»/T 




VhTft/2T 
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-fi/T 
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Following the argument used in deriving Eq. ( 13 .7j) . the entropy functional H 2 . 3 f) are generi- 
cally expanded as follows: 



AS(t; Z X [L S }) 



d D y 



[Ls 



e { ^) V^(ds/dp v) ) 
£(t d^ds/dn) 



1 ^2 (pa) 

2{nv), 



-3(pu),(pa) / 



y(p(ds/dp a) ) t 



pv gpv 1 

1 <-2 



1% If) \d v (ds/dn) 
(l\ t\ t\ I 



tre 9 V M (<9s/<%) 



Cs ps ps 
-2 l 3 c 5 



\ 



[y^ds/dp^J 



d D y 



(e { , v) {-l/T)K 



c #y W(-/vn 

( tr p \ 



+ (tre (-1/T)trif 



;7S /Js «s 

'2 c 3 c 5 

;7s *s ps 

-A c 5 c 6, 



^(-1/T)tr/^ 



(4.9) 
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where the coefficient matrices are symmetric and positive semidefinite, and their indices are 
taken to be all orthogonal to . Note that we again have used the fact that the rotation 
h p hj 7 V[pU CT ] vanishes here. Since the matrices must be invariant tensors, we can assume 
that they take the following form: 



^p,v),{pa) ^u),(po-y 



9 \ 4 \ 

<)i pi 

'2 l 3s 



$2 



'1 l 2 



(4.10) 
(4.11) 



where ' " 2 



md y t l q J are positive semidefinite. Note that 4^'^ 5 (p°-} = 
{k = 1, 2, 3) . The functional derivatives are then evaluated to be 

SAS 



5AS 
fen 
5AS 
6 (tie) 
5AS 



5AS' 
Se„. 



59 
5AS 



[l\ e<""> - (4/T) ] , 

^r-gh^[i\e v + ild u {-^/T)] 

S^g[l\ tre + f 2 9- (tjT) trK] , 

^[f 2 tTS + t 3 e-(£t/T) trK], 
N (d u*V + C || /i'*") x 

x V[4e M - (£l/T)K {up) + tr £ + £|0- (^/T) trif)/^] 



<5AS 



+ 



5 AS 



5AS 



-0 



V4^^^ + ^^^(-^/T)] 



where we have made the following parametrization as in Eq. (13.12j) : 

d 2 s 



dp^dp u 



- c± u»u v - en 



Note that SAS/de^ has only the components tangent to the timeslices. 
The equations of linear regression are given in the following form: 



1 re e SAS 



pa 



L Mbrr± ^ ± 



5AS 

5e v 



0, 



(4.12) 
(4.13) 
(4.14) 
(4.15) 



(4.16) 
(4.17) 

(4.18) 



(4.19) 
(4.20) 
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\e J . „ = -= Iff " 
L H irr || || 



5AS 



'[(t^)-] ir ; 

[Pl* - fe]rev]j 
[PM - [P/Jrev], 

n - fnl™ 



i / £tr e tr e £tr e t 

~7fi \ L Uee L ee 



f8AS/8(tre) 



[iVj irrl 
[^]irr|| 



'5AS 
. 

^AS 1 " 

Spu . 



n 



5AS 

8h 



(4.21) 
(4.22) 

(4.23) 
(4.24) 
(4.25) 



We here do not consider the direct coupling between the variables and , expecting that 
it is negligible. Note also that tie = = W v ' can be replaced by (tre)' at this order 

of approximation because (tre)' — tre = (h^)' e^ v = —INK^e^. As we have done for 
fluids, we make the following irreducible decompositions under the group 0(D) in a local 
inert ial frame: 

(4.26) 
(4.27) 
(4.28) 
(4.29) 
(4.30) 
(4.31) 
(4.32) 
(4.33) 

Thus, by assuming that the reversible evolution of and n is expressed with reversible 
currents rfi v and z/[\ (whose explicit form will be given later for a simple case) as 





= h'f h^hfii p h v 




= 17 hp V , 


t tr e tr e 


— T s 


^trefl 


— ^2 > 


L ee 


— T s 




= L± u^u u , 




= L\\ hp V , 




= M . 



'(r) 
%] 



n\ 



-9 V M ^ r) 



(4.34) 



Eqs. f l4^23|) -f l4T25|) become 



c± L± (-u u u x ) [4 6 { " x) - (1/T) 4 if^ A > 
+ (£\ tie + t 5 9- (1/T) £ s 6 tr K) h» x ] , 



(4.35) 



[V M (pV* + 7gU = - C |,L,| fc^V,, [4e<^> - (1/T) t\K^ 
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+ (£1 tre + f 5 e- {I /T) 1% tr K) h^] , (4.36) 
V^n^ + i/jy =-Vh(-8Ps/dn 2 )MV„[qh^e u + £ v 3 h^d u (- t i/T)]. (4.37) 

Assuming again that Lp = c±L± = c\\L\\ and L h = ^/h(—d 2 s/dn 2 )M are constant, the 
above equations can be written in the form of current conservations: 

V M T^ = , VX = , (4.38) 

where the conserved currents are given by 

T*» = e uV + 7$ + 7$ , n" = nu» + z/f r) + ^ , (4.39) 

with the dissipative currents 

= Lp 4 - (1/T) Lp 4 AT^ + L p - tr e + £\ 9 - (1/T) £ s 6 tr hT , (4.40) 
^ d) = L h [t 2 h^e v + qh^d v (-f,/T)] . (4.41) 

As for the b A -type variables, the equations of linear regression Eqs. (14 . 1 9 [) ( T4 . 2 2 p can be 
represented in the following form: 
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(4.45) 


£ u 9 
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L J rev 




V Llll)±txK. 








{L\ £\ + Lg 


3) tre- (L|4 + L s 3 4)0 + (L|£| - 


(4.46) 



They give the generally covariant extension of the rheology equations introduced in [i~0~l [TTj . 
and describe the dynamics of plastic deformations. 

The current conservations (I4.38l) -( l4.4ip and the rheology equations (I4.42| — (I4.46p totally 
determine the time evolution of the variables , n , , , and 9 , and thus constitute the 
set of the fundamental equations that govern the dynamics of relativistic viscoelastic mate- 
rials. Further studies of relativistic viscoelastic materials are performed in our subsequent 
paper [T5] . 
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We finally make a comment that the explicit forms of [£(^)] r ev > [^] r ' [(^ r£ ) ] > arL( ^ 
Lev depend on the system under consideration. For example, for the case when the 
reversible, isentropic evolution describes that of elastic materials, they are given by 

7$ = - 2G e W -K(ir£-a6) hT , = , (4.47) 
^" 1 [^>] w = ^>. ^[(treyl^^trK, (4.48) 
kJ„v=[*U = °- (4.49) 

Here the non-negative constants Q and K. represent the shear and bulk modulus, respec- 
tively, and the constant a is proportional to the coefficient of thermal expansion. Equations 
(14.481) and (14.491) show that the intrinsic metric does not vary for reversible processes 
( [^V] rcv = [ 2N ^^] rcv = 2AA' M „ - [e^] rcv = 0) . A few other examples of isentropic 
evolutions are also given in [T5j . 



5. Conclusion and discussions 



In this paper, we have developed an entropic formulation of relativistic continuum mechanics, 
which includes the standard relativistic fluid theory. The discussion is based on the linear 
nonequilibrium thermodynamics, and we have proposed a local functional S which represents 
the total entropy of nonequilibrium states and is to be maximized in the course of evolution. 
We have applied this framework to constructing a relativistic theory of viscoelastic materials. 
As is intensively studied in our subsequent paper [15], this theory can deal with a wide 
class of continuum materials, including as special cases elastic materials, Maxwell materials, 
Kelvin- Voigt materials, relativistic viscous fluids, and the so-called simplified Israel-Stewart 
fluids, and thus is expected to be the most universal description of single-component (not 
necessarily relativistic) continuum materials. 

Our entropy functional S is a local functional. Thus, once the coefficients £ AB , ■ ■ • , i" 1 ' vJ 
in Eq. (I2.30p are determined for a given material, one can explicitly evaluate the difference of 
the entropy of any configuration from that in the global equilibrium, AS — S — So , as long 
as the configuration is not far from the global equilibrium; one only needs to measure the 
local values of thermodynamic variables (such as the values of temperature at each point) 
and put the obtained data into AS . 

The next step would be to extend the current formalism as one can deal with more 
complicated systems like multicomponent viscoelastic materials or systems with extra vari- 



25 



ables like liquid crystals. Such extension is actually straightforward and is under current 
investigation. 

Another interesting direction would be to extract from our analysis the information on 
the holography of gravitation. In fact, one would need to assign local entropy to the metric 
tensor when the system is analyzed in generic frames (other than the Landau-Lifshitz frame) 
and when material particles have nonvanishing accelerations [16]. If this point is carefully 
investigated within the framework of the Bogoliubov-Born-Green-Kirkwood-Yvon hierarchy, 
we would be able to find out the microscopic degrees of freedom which need to be introduced 
to describe general relativity at large spacetime scales. Investigation along this line is now 
in progress, and will be reported elsewhere. 
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A. A derivation of the entropy functional 

In this appendix, we derive the entropy functional f)2.30p in the flat background, where 
the local entropy is assumed to take a simple form, s(x) = s(c 7 (x)j . In the following, we 
consider the special cases where A = 1 and hij are constant. 

As can be seen from the equilibrium condition f 1 2 . X 8 [) . all intensive parameters Pi are 
spatially constant in the equilibrium states. We take the region £ X [X S ] to be a D-dimensional 
torus with the coordinates x % having the period L s . We assume that finite size effects become 
irrelevant when L s ^> e s . Then the equilibrium values in H x [L s ] are equal to the mean values 
of c 7 (x) in T, X [L 8 ]: 

cl = ^- (V = VhL°), or % = (A.l) 
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where C 1 = j d D y c 1 (t,y) is the total charge for the region S Z [L B ] . We Fourier-expand 
c T (t,y) as 

~ cI (*> y) = v T, eiky %® = y + v^ eikv di{t) ' (A - 2) 

k fc^O 

where k{ = 2mii/ L s (71$ G Z) . The entropy functional 

S(t; S,[L S ]) = / 5(c 7 (t,y)) (A.3) 

«/ Sj; [Z/ S ] 

can then be expanded as follows: 

S = Vs(c I (t)) + ^ [ ^y^WYs^-clit) 

where s^j = (d 2 s/dc I dc J ) | Q , and we have used f d D y Xlfc^o elfc y ^fe(^) = ■ 
S — S can be written as 



7^ 

(v/i) 2 V 



2^, 



(A.5) 



with fu(k) = — s^j (1 — 5fc,o) • However, this nonanalytic form of fu(k) is not desirable 
because it may change the long-distance behavior in the analysis based on the derivative 
expansion. The most desirable is such a function that filters out low-fc modes analytically, 
and one can take as such the following function: 

fuik) I' ! ' j /.-, /.;, . (A.6) 

where £ i I J J is a positive definite tensor and is symmetric under the exchange (il) -H- (jJ) . 
The (mild) increase of fu(k) in the region |fe| > 2tt/L s does not cause a problem because 
c£(t) are assumed to decrease rapidly in the same region. Then by using the equality 
d^ds/dc 1 ) = (l/y/h) sijdid J , the entropy functional can be written as 

S = So ~ 2V ^ 2 £ tifj ** k i 

= J- / & D yt i ; j dic I (t,y)dfc J (t,y) 

27 



where £ lI ^ J = t K 3 L (sq ) KI (sq 1 ) lj . For a general metric, S will have the following form: 

s= s ° -5 Lf x NVKe "" J v "(§) v -(i) • (A 8) 

where the coefficients have only the spatial components, £i J - I > 1/J u u = . If there further exist 
additional b A -tjpe variables, this functional will have the more general form f)2.30p . 
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